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Green function Retrieval and Time-reversal in a Disordered World 
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We apply the theory of muhiple wave scattering to two contemporary, related topics: imaging 
with diffuse correlations and stability of time-reversal of diffuse waves, using equipartition, coherent 
backscattering and frequency speckles as fundamental concepts. 
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In its early days, multiple scattering of waves was con- 
sidered to be an unavoidable nuisance. It randomizes 
phase, polarization and wave vector of waves, and thus 
complicates important applications in imaging, telecom- 
munication, laser action, and remote sensing. Recent 
developments showed that multiple scattering can ac- 
tually enhance the performance of several applications. 
We mention in this context the low threshold of random 
lasers 0, the robust time-reversal of multiply scattered 
waves with potential applications in communication 
ISl, the sensitivity of diffuse waves to particle motions 
and the retrievals of the Green function from ther- 
mal phonons |^ and diffuse seismic waves ^i]. 

Chaos theory successfully describes time-reversal .7j 
and correlations ^ of elastic waves in chaotic media. 
Most applications above concern diffuse waves, with their 
specific questions about statistics, leakage and dynamics. 
It is the intention of this Letter to apply multiple scat- 
tering theory to these new exciting topics. To our knowl- 
edge the first attempts to cross-correlate "noisy" signals 
to retrieve ballistic wave motion come from seismology 
^acoustic daylight imaging") Q and helio-seismology 
As for time-reversal, after many pioneering exper- 
iments in Paris ^Jli Papanicolaou etal. explained 
the self-averaging property of time-reversal for broad- 
band signals for (locally) layered random media in a scal- 
ing limit of small fluctuations and propagation distances 
long compared to the wavelength. Our approach makes 
no assumptions on the scatterers, but applies only in the 
diffuse regime, i.e. when the medium size L is much big- 
ger than the mean free path I* . This excludes ID and 
localized media. This criterion is largely fulfilled in 2D 
time-reversal experiments with high order multiple scat- 
tering 0. 

We first consider the fiuctuations of a scalar wave 
field 5'(r,t) propagating in an infinite random medium 
after being released by some source S far away 
from the place of measurement. It is customary 
to describe field correlations by the "Wigner" func- 
tion (r- ix,i- ir) (r-l- ix, i-l- ir), whose 
Fourier transform with respect to space x and time r 
is defined as the specific intensity /c^k (r, t) of waves with 
wave vector k and frequency oj. A rigorous statement in 



multiple scattering theory is that 



lim /dr(/^k(r,t)) = --Im(G(^,k)) 
Eu'-^Im(G(c^,k'))5(c.,k') 



Ek'-^Im(G(c.,k')) 



(1) 



The structure function of the source S{Lo,'k) is the FT of 
5(t, x) = / dt / dr s (r - ix, t - \t) s{Y+\yi,t+ \t) , 
and (G(tJ,k)) is the average Green's function of the ran- 
dom medium. Equation is a known field-theoretical 
consequence of flux conservation 0, . It relates the 
speciflc intensity at large lapse times to the spectral func- 
tion — ^Im(G(w,k)) of the effective medium [l^. This 
implies global equipartition of average energy in phase 
space. 

Equation |Q looks like a manifestation of the 
fluctuation-dissipation (FD) theorem, applied to ^' = 
G ® s for which thermal equilibrium (TE) implies that 
{■^{u})-^*{uj)) oc Im {G{uj))B{u}), with B the Planck func- 
tion. In TE, formula ^ would even hold without mul- 
tiple scattering, and at any time. This situation applies 
to the experiment by Weaver and Lobkis Q, who mea- 
sured the elastic Green's function of an aluminium block 
by cross-correlating thermal phonons. When translated 
to space-time, the FD theorem shows that coherent wave 
paths should in principle be observable with travel times 
up to the inelastic phonon mean free time. 

The diffuse field is not in TE, not even in equilibrium. 
The physics that multiple scattering and TE share is the 
equipartition principle [3 ^| . Unfortunately, two rea- 
sons exist why Eq. ^ is not very useful to retrieve the 
Green function. First, it applies to the ensemble-average 
of the specific intensity only, which is subject to large 
statistical fiuctuations. This would clearly restrict its use 
in several imaging problems. Secondly, Eq. ^ assumes 
many detectors at different positions. As shown in Ref. 
5], the correlation method can be adapted for infinitely 
many equal sources, which also eliminates the problem 
of ensemble-averaging. In this work we will consider 
the random field generated by one distant source and 
investigate the applicability to retrieve the Green's func- 
tion using the local time correlation function <I>r (x, r) = 



2 



/ * (r - ix, i - ir) * (r + ix, t + ir) . We shall de- 
velop a theory for its statistical average, and show that 
<I>r(x,r) tends to be self-averaging. 

Equation suggests the following local expansion of 
the ensemble-averaged specific intensity, 



(/^k(r,i)) = --Im (G(c^,k)) 

UJ 



l~vX-k-dr 

UJ 



xS{Lu)pUr,t) .(2) 



where the energy density pi^(r,t) satisfies a diffusion 
equation with source d{r)6{t), and the power spectrum 
S{uj) is defined as the source factor on the righthand 
side of Eq. Q ■ The second term in Eq. ij^Jl gives a dif- 
fuse flow of energy and is known to be valid in dimen- 
sions d > 1. Let us consider an ideal explosive source 
s(r,T) = SS{t)S'{r), with power spectrum S{uj) — Slj^ 
that we band-filter in a frequency band B, inside which 
transport quantities like the diffusion constant D arc as- 
sumed not to vary too much. We get the simple result, 

($r(x, r)) = S [p{r)dr - vptdrpir) ■d^ + ---]x 

[(GB(x,r))-(GB(x,-T))]. (3) 

The ensemble-averaged field correlation function is pro- 
portional to the time-derivative of the band-filtered, 
ensemble- averaged Green's function (Gb) of the random 
medium, symmetrized in the time r [3- An antisymmet- 
ric part is allowed if a diffuse flux is present. There is a 
strong belief that this term is the origin of the asymmetry 
of the correlation function recently observed with seismic 
Rayleigh waves Q. The verification of this hypothesis 
requires a more realistic diffusion model for the Earth 
crust, including (free) surface detection, surface waves 
and mode conversions that is beyond the scope of 
this Letter, and for which laboratory experiments will be 
indispensable |23| . The diffuse character of these seismic 
waves had been deduced earlier from their equipartition 
PH and their seismic coda 22]. In an open medium of 
size L the relative importance of the asymmetric term 
is of order £* / L far away from its boundaries, and thus 
small if L > ^* . 

If the waves propagate in 3D without much dispersion, 
(G(x, r)) « —6{x — ct)/Attx for r small compared to the 
mean free time, and Eq. Q can easily be generalized for 
an arbitrary source spectrum. 
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S{t) is the time correlation function of the source. In par- 
ticular, the result ($r(0,T)) '-^ p{y)S{t) will be needed 
later to describe time-reversal. Relation Q facilitates 
the monitoring of source dynamics with distant, diffuse 
correlations. 

Equations 1)314(1 constitute the central result of this pa- 
per, despite their simplicity. Their importance follows 
from the rest of this Letter. First we establish that fluctu- 
ations around the ensemble-average are small if the band- 
width B is large enough, as also found for layered random 
media [T^ J l4 | . and first suggested by time- reversal ex- 
periments Q . We will generalize the central result Q to 
a random medium containing an arbitrary close object. 
This shows the usefulness of diffuse waves to the inverse 
problem in disordered media, in a way as first put for- 
ward by Claerbout etal. for seismic noise. At the end 
of this Letter we make the link with time-reversal and 
coherent backscattering in disordered media [TH . 

We start out by calculating the statistical fluctua- 
tions of $r(x,T) around its average. To this end let us 
adopt a source at a large distance r with power spec- 
trum S{u!) = Sbj"^ that we filter over a bandwidth B. 
Our most important assumption will be that the band- 
width B is much larger than the Thouless frequency 
rith — -D/2d|rp of the random medium (with D the dif- 
fusion constant, and d the dimension). It can readily be 
seen that (r='= = r ± ^x). 



Mx,r)=/ -e 



G(r-, 0, t^)G*(r+, 0, u)S{uj) .(5) 



in terms of the monochromatic, retarded Green's func- 
tion G(r, r',w) of the wave equation. From Eq. ||2Jl we 
find for the average correlation. 



($r(0,0)) 



Aoj -1 
27r UJ 



Im (G(w, G))Sp^{r) - BNSp{v) 



(4) 



with N the density of states in the bandwidth TS] . In the 
diffuse regime spatial /f req ucncy fluctuations obey Gaus- 
sian (Gi) statistics [23, so that the variance becomes. 



^ / ^ e-'("^-'^^)^5(c.i)^(c^2) (G(r-, 0, uJr)G* (r", 0, oj^)) (G(r+, 0, lo^) G*(r+, 0, oJ^)) 

Jb 27r 
~ B\S\''N^p{vfnt^. 
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We use that the normalized intensity correlation function 
C^(r2) ~ exp(— -^ri/ilth) decays rapidly with and that 
fith < B, 1/t. In that case is A$r/(^'r) - \/^th/B <C 
1, so that the measurement of <l>r(x, r) is close to its 
ensemble- average with high probability. The ratio B/ilth 
is interpreted as the number of independent "frequency 
bits" available in the bandwidth of the waves arriving at 
the receiver Q. 

In the following we will show the possibility to image 
a fixed, close object by means of diffuse correlations of 
a scalar field 5*. The formal relation between the Green 
function of the effective medium including the object, 
the Green function (G) of the effective medium alone, and 
the T-matrix T of the object is G'^ = G + GTG (we drop 
the averaging brackets if no confusion can arise). We 
assert that Eq. |j3Jl is universal so that at a large distance 
r from an explosive source, 

$r(x,r) oc/?(r) 

94(Gi(r-,r+,r))-(GS(r-,r+,-r))] (6) 

where B refers again to a finite bandwidth. The iden- 
tity © implies that the field correlation of the diffuse 
waves between two points near the object is equivalent 
to a time-resolved scattering experiment on the object. 

This assertion can be established from the following 
assumptions: 1. We are in the diffuse field of a distant 
(explosive) source S; 2. the object scatters the waves 
elastically. Our starting point is Eq. (0) valid without 
the object. The solution with object can formally be 
found by "gluing" a single scattering vertex K into the 
transport vertex p giving the terms {p+Kp+pK+pKp)S, 
with S the source. In momentum space the vertex K at 
frequency uj is given by |24| , 

+ ^k+iqi,k'+iq,G*,_\q^^k-iqi,k'-iq2 
-iqi,k'-iq2'^k'+iq/k+iqi,k'+iq2 > 



t: 



where we abbreviated G{uj,'k) — Gk. pKS is negligible 
for a distant source and pKpS just slightly modifies the 
diffuse background p. The terms pS + KpS dominate if 
the object is less than a mean free path separated from 
the receiver, and give us. 



*k+iq,*k-iq, 



qi^k-iqi 



V [5kk' + -ft:kk'(qi>q2)] — ImGk'/5(q2)5'(w) . (7) 



The Ward identity [2£ 



^P+^q,p' + iq ^p-iq,p'-iq ^ 



X!-^P+5q.k+iq 



^p'-iq,k-iq- 



1) imposes 52 ~ in all T-matrices and Green functions. 
The assertion © follows after Fourier transformations 
with respect to w,k, and q2. 

We will finally make the link of Eq. ^ with time- 
reversal (TR) experiments and coherent backscat- 
tering (CBS). Consider a source s(t) at point rg in an 
infinite random medium. At point a TR process is 
carried out between times T\ and T^. The TR signal 
^'(r, r) collected at ru was given by De Rosny etal. 
and Papanicolaou etal. (1^ : 

^'(rz5,r+2r2) = / dtG(rs,rT,T+i)G(rT,ri3,t)(^s(t,rs) 

(8) 

If we agree to time-reverse the whole signal (Ti smaller 
than the first arrival time and T2 deep inside the coda), 
a Fourier transformation gives us, 

/OO 1 
-00 27r 

xs(w) .(9) 

The similarity to <E'r(x, r) (with r = — rj- and x = 
rs — !■£)) given by Eq. (O expressing diffuse correlations 
is striking. Thus Eq. Q explains the auto-focalization 
in space-time, whereas the TR signal ^'(r, r) is very 
stable against statistical fluctuations provided it has a 
bandwidth much larger than the Thouless frequency. In 
Ref. the ratio equals B/fith « 600. 

The link between TR and CBS can be established by 
applying formula to a disordered half space with a 
TR machine hidden inside at a depth zt and a source at 
a large distance a outside (Figure 1). We assume zt -C a. 
Equation 1^ predicts that the TR signal at a distance x 
from the source is, 



«'(x,r + 2T2) = sir) y d^r,, e^'=''-'-|i/>(£, zt, r,,) 

ZTk\x 



s{t) — exp 

ZT 



(10) 



where we have inserted the diffuse e nerg y propagator 



expressing flux conservation (assumption 2) simplifies 
this expression. The singular diffusion pole (assumption 



p(z,z', r|| — rj|) for the 3D half space |2g. This makes 
the line shape of the spatial auto-focalization in TR equal 
to the one of CBS from a half space, exhibiting the 
well-known angular triangular cusp caused by very long 
time-reversed wave paths 26]. The angular resolution 
Ax /a « X/zt nicely illustrates that the diffusion pro- 
cess creates an effective aperture of size zt that increases 
the quality of spatial focalization, independent of i* 0]. 
The peak signal is stable against statistical fluctuations 
if B » 2D/z^. 

In conclusion, we have shown that field-correlations of 
diffuse waves can be used to retrieve ballistic waves be- 
tween two points in space-time. This method is stable 
against the unknown statistical fluctuations that cause 
the field to be diffuse, provided that the operating band- 
width is much larger than the Thouless frequency. We 
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FIG. 1: A set-up that establishes the Unk between coherent 
backscattcring and time-reversaL A source S emits a signal 
that is received by a time-reversal machine T inside a disor- 
dered half-space. The signal that is sent back auto-focalizes in 
space and time near S with the spatial line-shape of coherent 
backscattering. 



have discussed the possible temporal asymmetry of the 
field correlations, the role of the power spectrum of the 
source, as well as the fundamental relation to time- 
reversal and coherent backscattering. The method facil- 
itates a passive way of imaging that might find an appli- 
cation in seismology, where active sources are expensive 
or unpredictable. 
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